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Abstract 

We perform linear stability analyses (LSA) and direct numerical simulations (DNS) to investi¬ 
gate the influence of the dynamic viscosity on viscous fingering (VF) instability in miscible slices. 
Selecting the characteristic scales appropriately the importance of the magnitude of the dynamic 
viscosity of individual fluids on VF in miscible slice has been shown in the context of the transient 
interfacial tension. Further, we have confirmed this result for immiscible fluids and manifest the 
similarities between VF in immiscible and miscible slices with transient interfacial tension. In a 
more general setting, the findings of this pletter will be very useful for multiphase viscous flow, in 
which the momentum balance equation contains an additional stress term free from the dynamic 
viscosity. 
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Displacement processes through porous rocks and mixing of two miscible fluids are active 
areas of research, having several industrial and environmental applications, such as enhanced 
oil recovery, hydrology and filtration, carbon capture and storage, etc. VF, a hydrodynamic 
instability that occurs in both the immiscible and miscible fluids while displacing a more 
viscous fluid by a less viscous one, is inherent in such flow configuration DH3j. In immiscible 
fluids surface tension force at the interface acts against the instability [IjJ. On the other hand, 
in miscible fluids, where a thermodynamically stable interface does not exist, a transition 
zone relaxes with time due to diffusion and acts against the finger growth. Experiments EE] 
reveal, when the diffusion is slow, a steep gradient in the form of density, concentration or 
temperature between the underlying fluids gives rise to a weak transient interfacial tension 
that mimic surface tension effect. This was first discussed by Korteweg in 1901 [Tj, who 
introduced an additional stress term, known as the Korteweg stress, in the equation of 
motion. The existence of Korteweg stress or transient surface tension is also observed in the 
experiments of colloidal suspensions 0 and in the binary liquid system of isobuteric acid 
and water [9j. 

Chen and Wang analyzed the influence of VF instability on the spreading of a local¬ 
ized fluid slice having higher mobility than the surrounding fluid. On the other hand, De 
Wit et al. mi studied the same problem in the context of separation in a chromatographic 
column when the viscosity of the sample is higher than the solvent. Mishra et al. ca have 
shown that the onset of VF instability and the subsequent finger pattern near the onset are 
identical for both the less and more viscous slice. Influence of the Korteweg stresses on VF 
instability was investigated theoretically by Joseph and his co-workers [131 iQJ and Chen 
et al. na However, to the best of the authors’ knowledge the influence of such stresses 
on the nonlinear VF instability of more and less viscous miscible slices in a Helc-Shaw cell 
has never been addressed adequately. In particular this letter addresses the question, what 
is the influence of the Korteweg stress that describes volume forces arising because of the 
nonlocal molecular interactions on the VF at the rear and frontal interfaces of a localized 
slice? Such a classical complex pattern dynamics has been investigated through a highly ac¬ 
curate Fourier-spectral method based direct numerical simulations m- It has been proved 
theoretically that an appropriate choice of the dynamic viscosity of underlying fluids results 
the identical onset of fingering instability and the subsequent finger patterns are also iden¬ 
tical for both more and less viscous slices in the presence of the Korteweg stresses. The 
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dc / By = 0, dv / dy = 0 
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dc / dy = 0, dv / dy = 0 

FIG. 1: Schematic of the flow in a Hele-Shaw cell. 


DNS results are found to be in excellent agreement with the corresponding LSA. Also, the 
similarities between the immiscible slices and miscible one with transient interfacial tension 
have been proved through LSA, which affirms the classical nature of this study. 


I. MATHEMATICAL FORMULATION AND NONLINEAR SIMULATIONS 

Consider a uniform rectilinear displacement of a finite fluid slice of viscosity p 2 by another 
fluid of viscosity p i in 2D porous media or a Hele-Shaw cell as shown in fig. [l] The frontal 
interface becomes unstable if pi < p 2 , otherwise it is the rear interface which features the 
fingering instability, for pi > p 2 . The viscosity of the fluids depends on the solute concen¬ 
tration c, i.e. p = p(c). Fluids are assumed to be incompressible and neutrally buoyant. 
With an additional condition of slow diffusion, the above-mentioned flow problem can be 
described in terms of Darcy-Korteweg equations USE] coupled with a convection-diffusion 
equation for the mass conservation of the solute concentration. For the dimensionless formu¬ 
lation of the equations the diffusive length and time scales, D/U and D/U 2 , are used as the 
respective characteristic scales. Characteristic pressure, velocity, concentration and viscosity 
are taken to be p\ D/k, U, c 2 and p i, respectively. Here n is the constant permeability of the 
homogeneous porous medium. For simplicity we have assumed a constant isotropic diffusion 
of the solute concentration, characterized by the diffusion coefficient D. The dimensionless 
equations in a Lagrangian frame of reference, moving with the speed U are written as, 

V • u — 0, Vp — —/i(c)(u + i) + 5 V • (Vc <g) Vc), (1) 

c t + u ■ Vc = V 2 c. (2) 

Here u is the gap-averaged velocity having longitudinal and transverse components u and v, 
respectively, p is the dynamic pressure, 6 = 5nU 2 c\/piD 3 (< 0) is the dimensionless Korteweg 
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stress constant US ns and the operator V = iJ^ + The governing equations dlh - 
Q are associated with the following boundary conditions: at the longitudinal boundaries, 
u = (0,0 ),dc/dx = 0,x — > ±oo; and at the transverse boundaries, dv/dy = 0 (representing 
the constant pressure), dc/dy = 0,Va;, in the Lagrangian frame of reference (shown in 
fig. g. The initial velocity is considered to be u = (0,0), while the initial distribution 
of the solute concentration is, c = 1 inside the finite slice and c = 0 outside of that. 
The relationship between the dynamic viscosity of the underlying fluids and the solute 
concentration is assumed to be of Arrhenius type, y(c) = e Rc , where R = ln(/i 2 /jUi) is the 
log-mobility ratio mini. Hence, the displacement of a less (more) viscous slice by a more 
(less) viscous ambient fluid is represented by R < 0 (R > 0). 

Direct numerical simulations of eqs. ( 0-0 in terms of the stream function, ij}(x,y,t) 
(u = il) y ,v = -^x), 
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V 2 ^ = —RVc ■ (V^ + j ) H—y- Vc x V (V 2 c) • k, 


h(c) L 


Ct T 'IpyCx 'IpxCy 


( 3 ) 

( 4 ) 


are performed using a highly accurate pseudo-spectral method mi. Here k represents a 
unit vector normal to the xy-plane. Numerical simulations are performed in a rectangular 
domain of dimensionless width L y = 1024, and length A ■ L y = 8192, where A = L/H , 
the ratio of the length to the width of the domain, is the aspect ratio. A convergence 
analysis of the numerical method and grid independence have been performed. Nonlinear 
simulations are carried out in the above-mentioned computational domain with 2048 x 128 
and 2048 x 256 spectra, and it has been observed that the maximum relative error between 
these two sets of simulations is of O(10~ 3 ). Further refinement of the spatial disretization 
points in either direction does not alter the dynamics of the fingering pattern. Hence, for 
optimal computational cost we discretize our computation domain of size [0, 8192] x [0,1024] 
using 2048 x 128 spectral points. 


II. RESULTS AND DISCUSSION 

In the absence of the Korteweg stress VF dynamics for both R > 0 and R < 0 are 
the same until the interaction between the two interfaces [T2] , Our aim in this letter is to 
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(a) (b) 



FIG. 2: Density plots of concentration at successive times in the Lagrangian reference 
frame for L y = 1024, l = 256 ,5 = —10 4 when (a) /i = e 3c , (b) /i = e _3c . From top to 

bottom t = 500,1000,1500. 


determine the possible changes in the dynamics of such slices in the presence of the Korteweg 
stress and relate it to the surface tension effect in immiscible slices. DNS results for \R\ — 3, 
L y = 1024, Z = 256, <5 = —10 4 are shown in fig. @ which depicts that at t — 500 fingers 
appear at the rear interface of the slice for R = 3 (see fig. [2j(a)) , whereas, for R = — 3 the 
dynamics of both the interfaces feature diffusive characteristic, but no finger formation at 
the frontal interface having the unstable viscosity gradient (fig. [2](b)) . This signifies that 
the presence of the Korteweg stresses induces different onset of finger formation for R > 0 
and R < 0. Moreover, the wavelength of the unstable modes is larger for R < 0 than that 
corresponding to R > 0. These differences between R > 0 and R < 0 persist in the long-time 
behavior of the fingering dynamics of the slices. 

It is customary to ask what causes such differences in the dynamics in the presence of 
the Korteweg stress? Conventionally, the fi — c relation, /j = e Rc , keeps ^ 2/^1 unchanged 
for both more and less viscous finite slices, but the dimensionless value of the less (more) 
viscous fluid changes with the sign of R. We notice that in the absence of transient interfacial 
tension eq. ([3]) is free from /x(c) (see eq. (9) in [12] )• Hence, the onset of fingering and the 
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FIG. 3: Density plots of concentration at successive times, t = 500 (top), t = 1000 
(bottom), in the Lagrangian reference frame for L y = 1024,/ = 256: (a) 

5 = —5 x 10 2 ,p, = e~ 3c , (b) 6 = -10 4 ,/x = e 3 M. 


qualitative as well as the quantitative pattern of the fingers are identical for R > 0 and 
R < 0 in the absence of the Korteweg stress. 

Choosing the displacing fluid viscosity, /j j, as the characteristic scale for the dynamic 
viscosity in both the cases of more and less viscous slices we should be more careful about 
the dimensionless Korteweg stress constant, S , which is inversely proportional to the char¬ 
acteristic viscosity. As an example, consider the displacement of a glycerin-water mixture 
by pure water. The dynamic viscosity of glycerin-water mixture with 68% glycerin con¬ 
centration is 19.4 cP (~ 20 cP) and that of pure water is 1.002 cP (~ 1 cP). Therefore, 
in this case we have /xi ~ 1 cP and /r 2 ~ 20 cP, which can be written in terms of the 
dimensionless parameter R = 3. For this pair of miscible viscous fluids, a dimensional value 
6 = —10 -10 N corresponds to the dimensionless value 5 = —10 4 for an injection speed 
U ~ 10 -3 ms -1 and diffusivity D = 10 -9 m 2 s _1 . On the other hand, for the displacement 
of water by glycerin-water mixture ( C g = 68%) /i\ ~ 20 cP and /j ,2 ~ 1 cP. This makes the 
dimensionless value 5 = —5 x 10 2 for the same dimensional value 5 = —10” 10 N with the 
same characteristic variables. In order to justify our claim we perform numerical simulations 
for R = —3,5 = —5 x 10 2 , and the results obtained are presented graphically in fig. [3](a) , 
which depicts that the onset of instability and the finger pattern, until the interaction of 
the unstable and the stable interfaces, are identical to those for R — 3, 6 — —10 4 (fig. [2](a) ). 
This suggests, in order to compare the numerical results of two similar flow configurations 
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one must choose the dimensionless parameters in such a way that they correspond to the 
same dimensional value. Subsequently, we ask the question whether there exists a suitable 
characteristic scale for the dynamic viscosity which automatically takes care of this fact. 
In particular we choose the less viscosity /p as the characteristic viscosity and \xiD/k as 


the characteristic pressure (see table ??). Hence, the dimensionless form of fi — c relation 
becomes, 


H(c) = e R ’^\ /(c) = 


C, 


hi < h2 


R = 


1 - C, /lx > h2, 

In (h2/hi), hi < h 2 
hi (hl/h2) , hi > h2- 


(5) 

( 6 ) 


With this modifications we have /? > 0 for both the more and less viscous slices, and the 
initial dimensionless dynamic viscosity of the ambient fluid (finite slice) becomes 1 and e R 


( e R and 1), respectively (see table ??). Thus the magnitude of the characteristic viscosity, 


hchar’ rema i ns the same for both the less and more viscous slices that results the same S 
corresponding to every single dimensional value 5, in both the cases. Since, the mathematical 
description of the displacement of a high viscous slice by a low viscous ambient fluid remains 
the same after the modification, in what follows is the discussion of the displacement of a 
low viscous slice by a high viscous ambient fluid to understand the consequences of such 
modifications on the VF instability of more and less viscous slices in the presence of the 
Korteweg stresses. Fig. DM shows the concentration distributions at successive times for 
// = e 3 ( 1-c \ L y = 1024, Z = 256 ,5 = —10 4 , which are identical to those in fig. [3|a), hence 
depicts that the onset of VF and finger pattern are the same to those for a more viscous 
slice (fig. |2](a)) , until the interaction of the fingers with the respective stable interface. 

The onset of fingering instability and the interaction between the stable and unstable 
interfaces can be quantified appropriately from the temporal evolution of the interfacial 
length [Oi /(f) = f 0 L »f 0 ALy (it) 2 + (jt) dxdy- In the diffusion dominated regime 
/(f) retains a constant value equal to the width of the domain, and increases with the 
growth of the fingers. Hence, the onset of fingering is marked as the instance when /(f) 
starts increasing from the constant value of the diffusive regime. Fig. [4] shows the temporal 
evolution of the interfacial length corresponding to the simulations of figs. [2] and [3j It 
depicts that the displacements remain stable over the first two decades of time and the 
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Ml(cP) /« 2 (cP) 

d'char 

hf 

hf 

R 

1.002 

19.4 

hi 

1 

19.3613 

ps 3 

19.4 

1.002 

hi 

1/19.3613 

1 ! 

ps -3 

19.4 

1.002 

h2 

19.3613 

1 

ps 3 


TABLE I: Dimensional and dimensionless (superscripted ‘d’) values of the dynamics 
viscosity of the underlying fluids, characteristic viscosity, p ( .} iar , and R used in the 

computations. 



10* 10% io 3 

FIG. 4: Interfacial length I(t) for // = e 3c , (-) /i = e~ 3c (-) and // = e 3(1_c ) 

(-). The inset images show the concentration distribution approximately at time of 

interactions between the stable and unstable interfaces. 

fingers form at t ~ 180, for the viscosity relations, //(c) = e 3c and //(c) = e 3 *- 1-c y while at 
t ps 800 for //(c ) = e~ 3c . I(t) decreases with time after the interaction of the fingers with the 
respective stable interface. Inset images show the concentration distributions approximately 
at such instance of interaction and depict how this time differ from each other for different 
viscosity profiles. Hence, all the qualitative features of the dynamics observed in the absence 
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of any gradient stress [12], can be seen even with the gradient stresses, but with different 
quantitative measures due to the influence of transient interfacial tension. 

(a) (b) 




FIG. 5: Temporal evolution of y(t) for a more and less viscous finite slice of width, l = 256 

for L y = 1024 and 8. 

VF enhances the concentration gradient and the area available for diffusive flux across 
the interface by stretching the fluid-fluid interface and hence modifies the rate at which 
mixing occurs. It would be interesting to understand the effect of the Korteweg stresses 
on the mixing of a more or less viscous finite slice. The degree of mixing is defined as 
[T9] , y(t) = 1 — a 2 (f)/af liax , in terms of the global variance of the concentration field, 
cr 2 = (c 2 ) — (c) 2 , where (•) represents the spatial averaging over the domain. The degree 
of mixing of more and less viscous finite slices have been plotted in fig. [5] corresponding to 
the simulations of figs. 1-1 Due to stronger instability for the case of R — 3,8 — — 10 4 
than R = —3,8 = —10 4 , the degree of mixing, y(t), is higher in the former than the latter 
(see fig. Ha)). In fig. gb), x(t) has been shown for the three cases: (i) fi = e 3c ,8 = —10 4 , 
(ii) fi = e~ 3c ,8 = —5 x 10 2 and (iii) /i = e 3iyl ^ c \ 8 = —10 4 . It depicts identical temporal 
evolution of \{t) for both more and less viscous slices until the interaction of the fingers 
with the respective stable interface. Afterwards, y(f) becomes higher in the latter case than 
the former, and it asymptotically saturates to the maximum value of 1, corresponding to 
the completely mixed state. The coincidence of the three curves in the time interval t ~ 650 
confirms the results shown in figs. HN.i 
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III. LINEAR STABILITY ANALYSIS 


Here we discuss LSA of miscible slice in the presence of the Korteweg stresses with the 
viscosity relations given by eqs. ([5])-([6]) and compare the results obtained with those of DNS. 
Finally, the obtained LSA results are compared with a LSA of immiscible slice which confirms 
that the Korteweg stress and the surface tension have identical effects on the instability. 


(a) (b) 



FIG. 6: (a) Dispersion curves at to = 20, (b) Temporal evolution of the growth rate of the 
most dangerous wave perturbation (k m « 0.15). Inset image shows the maximum growth, 
amax^o, from all possible wave number, k, representing the onset of instability. The 
parameter values are 5 = —10 3 ,i? = 3 and l = 100. 


A. Miscible slice 

The initial-boundary value problem described above possesses a self-similar diffusive 
decaying solution in a similarity transformation (£, £)-domain [16], Uo = (0,0), /io = 
H o(£), Co(£) = [erf(£/2) — erf((£ — /*)/2)]/2, where £ = x/y/t is the similarity variable. 
Here, C*o(£) corresponds to the self-similar decay of a rectangular function of width l* = l/y/t 
and height unity. We analyze how an infinitesimal perturbation to this reference solution 
grows (decays) with time. Eliminating the pressure and transverse velocity component from 
the linearized form of the eqs. and using the fi — c relation from eqs. (§-(§, a 

system of coupled partial differential equations is obtained. These equations are homoge¬ 
neous in y and t, hence we can decompose the perturbation quantities into normal modes 
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(a) (b) 



FIG. 7: (a) Temporal evolution of the growth rate for 

R = 3,l = 100, k = 0.15, n = e Rc (0),fi = e R M (□), 8 = -10 3 . IVC (-) and DNS 

(-). Inset image is magnified near the onset of instability, (b) Dispersion curves at 

t 0 = 20 for l = 100 with the self-similar LSA. 

as, (v! ,d)(£,y,t) = ($(£), vi/(^))e* fcy+CT * (fc ’ to ^. Here T(£) and $(£) are the amplitude of the 
perturbed concentration c' and axial velocity v!, respectively, with k, a*{k, t 0 ) being the wave 
number and the growth rate of the perturbations [16] . Hence, the linear stability problem 
can be written as the following eigenvalue problem, 

(d 2 + f'(c 0 )DC 0 D - k%y = (^f'(c 0 )k 2 t 0 

+(k 2 S/y 0 Vi~o){D 3 C 0 - DC 0 (D 2 - k \) })*, (7) 

- (l/t 0 )^ 2 + k 2 - (£/2t 0 )D)V = -DCo®, (8) 

where D n = d n /d^ n ,?z G N. The system of coupled ordinary differential equations (eqs. 
([7])-(j8])) has been solved using a finite difference method [16] to determine the instantaneous 
growth rate of the perturbations in terms of the wave number k and the frozen diffusive time 
to- Dispersion curves obtained from LSA for the displacement of more and less viscous fluid 
slice are presented in fig. [6^a) for the parameters R = 3,1 = 100, <5 = —10 3 at the frozen 
diffusive time f 0 = 20. It depicts that the two dispersion curves are almost indistinguishable 
for the unstable modes. Also, the temporal evolution of the growth rate of the most unstable 
mode, k rn = 0.15 (fig. |6](a)), as well as the onset of instability for all possible modes coincide 
for n = e 3c and y = e 3 *- 1- ^ (fig. [6](b) ). Linearizing the stream-function form of eqs. 0-@ 
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with ii — c relation given by eqs. (|5])-(j6]) and using a pseudo-spectral method, an initial 
value calculation (IVC) is perfromed. The growth rate of the perturbed concentration, 
d, a c j = (l/2E c i)(dE c i/dt), where E c > = f^ v Ly d 2 (x,y,t)dxdy, obtained from IVC is 
compared with those obtained from DNS EDI. The temporal evolution of the growth rates 
of IVC coincide with those obtained from DNS for both /i = e 3c and // = e 3 ^ 1 ”^ with 
5 = —10 3 (fig. |r](a)). The numerical algorithm of the IVC will be discussed elsewhere in 
detail. Fig. [7](b) depicts an excellent agreement between the dispersion curves obtained from 
self-similar LSA at t 0 = 20 corresponding to R — 3, 8 — —10 4 and R = —3,5 = — 5 x 10 2 , 
and confirms the observation from DNS presented in fig. |3](a). 



FIG. 8: Dispersion curves for a three layer immiscible Hele-Shaw flow with 
L — 5, S — 1, T — 1, U = 1 for three sets of viscosities: /i/ = y r = e 2 , y = 1 (-), 

Ht = Hr = 1, n = e“ 2 (-), Hi = Hr = 1, n = e 2 (Q)- 
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B. Immiscible slice 


In what follows a comparison between an immiscible slice and a miscible slice with the 
Korteweg stresses. Chen et al. [ 21 ] have shown that in a rotating Hclc-Shaw cell both the 
qualitative and quantitative features of the Korteweg stress in miscible fluids and surface 
tension in immiscible fluids are identical. Not only the stabilizing property of the Korteweg 
stresses was confirmed, but also the number of fingers were measured to be the same for the 
miscible and immiscible fluids, both theoretically as well as experimentally. Our aim is to 
understand the relative importance of the viscous force to the surface tension force in three 
layer immiscible displacement, in which the growth rates of the perturbations at the two 
interfaces can be represented as (see eq. (10) in [ 22 ]), 0± = (—6 =F Vb 2 — 4ac)/2a, where, 
a = -e fci ( y u+ / u z )(/i+/i r .)+e _fci (/r- / u/)( / u-/i r .) < 0, b = [e fcL (/x+/x;)+e _fcZ '(/x-/r;)]^-[e fcZ '(/r+ 
/i r )+e _fcL (/i-/i r )]r, c = r£,(e kL -e~ kL ), £ = Tk((fi r -/j,)U/T-k 2 ), r = Sk(k 2 -(fi-fii)U/S). 
Here cr + and cr_ correspond to the larger and smaller growth rates, respectively. Furthermore, 
S and T are the surface tension at the left and right interfaces, L being the width of the 
middle layer, U is the fluid velocity and /q,/i, fi r are the viscosities of the left, middle and 
right fluid layers, respectively. These show that the obtained analytic expressions for the 
dispersion relations at the two interfaces involve explicitly the dynamic viscosities of the 
three fluid layers, not their ratio. As a consequence, the dispersion curves corresponding to 
cr_|_ become different for two sets of viscosities although having the same viscosity ratio of 
two adjacent fluid layers. For instance, we consider three sets of viscosities in such a way 
that the interface separating the middle layer from the right layer becomes unstable for the 
first two cases, while for the third, it is the interface between the left and the middle layer 
that features instability. Corresponding dispersion curves (cr + vs. k) are shown in fig. [8] for 
L — 5, S — 1, T — 1, U — 1; the solid line corresponds to the case of /q = /i r = e 2 , n = 1, 
dashed line corresponding to Hi — Hr — 1, /i — e~ 2 , and the circles representing the situation 
when hi = /i r = 1, h = e 2 . The viscosities of the left and right layers are chosen to be equal so 
that it reduces to an analogous form of a finite miscible slice discussed above. Although the 
ratio of the high viscosity to the low viscosity remain to be the same, e 2 , the dispersion curve 
for the second case differs from the remaining two. It signifies that the dispersion curves 
corresponding to a more and less viscous middle layer will overlap if the dynamic viscosities 
of the less and more viscous fluids are kept unchanged. Eventually, growth rate is always 
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higher for the case with lager dynamic viscosity of the fluids. Another important observation 
is the variation of the cutoff wave number (k c ) for the unstable mode and the most unstable 
wave number (k m ) with the variation of the surface tension. The qualitative changes of these 
two quantities are similar to those observed in miscible displacements under the influence 
of the Korteweg stress. Both in miscible and immiscible fluid systems, k c and k m decrease 
with the increase of the magnitude of S and the surface tension force, respectively. 


IV. CONCLUSION 


We have theoretically investigated the VF instability of miscible slices with transient in¬ 
terfacial tension. An appropriate choice of the dynamic viscosity of the fluid in the charac¬ 
teristic scales manifests its importance in the study of VF with transient interfacial tension. 
Our analysis is capable of reproducing the results existing in the literature. The results 
obtained depict that the onset of VF and their dynamical pattern can become identical only 
under certain scaling analysis for both the more and less viscosity slices in the presence of 
Korteweg stress, similar to the case when these stresses are absent ra- Synergetic mixing 
with VF and alternating injection |23j will be the same with the present viscosity model 
irrespective of the choice of the viscosity of the fluid that fills the Hcle-Shaw cell before the 
injection starts. Alike effect of the dynamic viscosity is also observed in immiscible fluids. 
The findings of this letter will certainly help to understand multiphase viscous flows with 
different viscosity of each phase or when the viscosity depends non-monotonically on the 
solute concentration |24j : for instance, during the mixing process of chemical components, 
pollutant contamination in aquifers, CO 2 sequestration, etc. Miscible displacement of vis¬ 
cous fluids with Korteweg stress has important applications in chemistry m- In this letter 
we present a preliminary understanding of the relative importance of the Korteweg stress 
to the viscous stress by incorporating one of the simplest situation of constant diffusivity, 
D, of the solute concentration. A velocity induced or concentration dependent diffusivity 
should be studied for a more realistic situation. More interesting and challenging problem 
of miscible VF in the presence of the Korteweg stress induced by both the concentration 
and temperature gradients is the focus of the authors’ future work. 
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